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The Dirac equation is generalized io D + \ space-time. The conserved 
angular momentum operators and their quantum numbers are discussed. 
The eigenfunctions of the total angular momenta are calculated for both 
odd D and even D cases. The radial equations for a spherically symmetric 
system are derived. The exact solutions for the system with a Coulomb 
potential are obtained analytically. The energy levels and the corresponding 
fine structure are also presented. 
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I. INTRODUCTION 



The exact solutions of the Schrodinger equation in the real three-dimensional space 
for a hydrogen atom and for a harmonic oscillator were important technical achievements 
in quantum mechanics During the past half century, the mathematical tools for 

the orbital angular momentum operators and their eigenfunctions in an arbitrary D- 
dimensional space have been presented . The nonrelativistic D-dimensional 

Coulombic and the harmonic oscillator problems have been studied in some detail by 
many authors [^,|lO|,|lT],|l2|,|l3|,|l4|,|l5l. The solutions of the Dirac equation, however, have 
been studied in the usual three- |ll7| , p!8| , p!9| , |20| , pT| , two- p2[ and one-dimensional p3| space. 
Motivated by the recent interest of higher- dimensional field theory , we generalize the 
Dirac equation to D+1 space-time. The conserved total angular momentum operators and 
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their quantum numbers are discussed. The eigenfunctions of the total angular momenta 
are calculated for both odd D and even D cases. From the viewpoint of mathematics, this 
problem is a typical application of group theory to physics. In terms of the eigenfunctions, 
we obtain the radial equations for the spherically symmetric system, and analytically solve 
the radial equations for the quantum Coulombic system. 

This paper is organized as follows. Section 2 is devoted to the generalization of the 
Dirac equation to D + 1 space-time. In Sec. 3, the conserved angular momentum oper- 
ators and their quantum numbers are discussed. The eigenfunctions of the total angular 
momentums are calculated for both odd D and even D cases in terms of the method of 
group theory. The radial equations for the system with a spherically symmetric potential 
are derived. In Sec. 4, the wave functions of bound states for a Coulombic system, which 
are expressed by the confluent hypergeometric functions, are presented. The energy levels 
and the corresponding fine structure are also discussed. Some conclusions are given in 
Sec. 5. 

II. DIRAC EQUATION IN I? + 1 DIMENSIONS 

The Dirac equation m. D + 1 dimensions can be expressed as 



At=0 

where M is the mass of the particle, and (D+1) matrices 7^ satisfy the anti- 
commutative relations: 

T^T^' + tV = Sr/^^'l, (2) 
with the metric tensor r/^'^ satisfying 

= = I "^/^^ ^^^^ = (3) 
\ -Sfj,u when /i 7^ 0. 

For simplicity, the natural units h = c = 1 are employed throughout this paper, 
consider the special case where only the zero component of is non-vanishing and 
spherically symmetric: 

eAo = V{r), Aa = 0, when a 0, (4) 
The Hamiltonian H{x.) of the system is expressed as 

D 

z9o^(x,t) = if(x)^(x,t), H{^) = J2 7VPa + l^(r)+7°M, 

Q »=i (5) 

Pa = -ida = -iTT-z, 1 <a< D. 

OX°' 
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The orbital angular momentum operators Lab, the spinor operators Sab, and the total 
angular momentum operators Jab are defined as follows 

Lab = -Lba = iXadb - iXbOa, Sab = -Sba = I'^alb/'^, 

Jab = Lab + Sab, 1 < a < b < D , ■ 

D D D (6) 

J^ = ^ Jab-: = X! -^afe' "S*^ = ^ 5*^^. 

a<b=2 a<b=2 a<b=2 

The eigenvalue of J^ (L^ or S*^) is denoted by the Casimir C2(M), where M is the 
highest weight of the representation to which the total (orbital or spinor) wave function 
belongs. We will discuss the Casimir in the next section. It is easy to show by the 



standard method that Jab and k are commutant with the Hamiltonian -ff(x), 



„o 



= 7" I E ^7 V^a. + {D-l)/2\=^'{j'-L'-S' + {D-l)/2}. (7) 

[a<b J 



III. THE RADIAL EQUATIONS 

Since the potential V{r) is spherically symmetric, the symmetry group of the system 
is SO{D) group. Erdelyi 0] and Louck introduced the hyperspherical coordinates in 
the real £)-dimensional space 

= r cos^^i sin ^2 • • • sin^^/)_i, 
= r sin 6i sin 62 ■ ■ ■ sin 6]j_i, 

x^ = r cos 9b-i sin 9k ■ ■ ■ sin 9d-i, 3 < b < D — 1, 

Do ^ ' 

X = r cos 9d^i 

a=l 

The unit vector along x is usually denoted by x = x/r. The volume element of the 
configuration space is 

D D-1 

II dx" = r^'^drdn, dn= H (sin d9a, 

a=l a=l 

< r < 00, -n <9i<Tr, 0<9b<n, 2 <b < D - I. 

Now, let us sketch some necessary knowledge of the SO(D) group. From the represen- 
tation theory of Lie groups P5|PB| , P7| , the Lie algebras of the S0(2iV+l) group and the 
S0(2A^) group are Bn and D^, respectively. Their Chevalley bases with the subscript /i, 
1 < /i < iV — 1, are same: 

H^{J) = J{2fi-l)(2fj.) — Ji2fi+l)(2fj.+2), 

Ef^iJ) = ('/(2m)(2^+1) - iJ{2^i^l){2ti+l) - 'iJ(2fi){2ti+2) - </(2/i-l){2M+2) ) /2, (lOo) 
Ff,{J) = (^(2/i){2M+l) + iJ{2f,^l){2,i+l) + iJ{2,i){2f,+2) " J{2f,-l){2,i+2)) /2. 
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But, the bases with the subscript are different: 

Hn{J) = 2J(2Ar-l)(2Af), 

En{J) = J{2N){2N+l) — iJ{2N-l){2N+l), 

(106) 

Fn{J) = J(2N){2N+1) + 'iJ{2N-l){2N+l), 

for S0(2iV + 1), and 

Hn{J) = J{2N-'i){2N-2) + J{2N-1){2N), 

En{J) = (^J{2N~2){2N-1) — iJ{2N-3){2N-l) + iJ(2N-2)(2N) + J(2N-i){2N)) /2, (lOc) 
Fn{J) = (yJ {2N -2){2N -I) + 'iJ{2N-?,){2N-l) — {2N ~2){2N) + J (2N -•i)(2N)j I'^i 

for S0(2A^). The operator Jab can be replaced with Lab or Sab depending on the wave 
functions one is discussing. H^{J) span the Cartan subalgebra, and their eigenvalues for 
an eigenstate |m) in a given irreducible representation are the components of a weight 
vector m = (mi, . . . , rriTv): 

iJ^(J)|m) =m^|m), l<^l<N. (11) 

If the eigenstates |m) for a given weight m are degenerate, this weight is called a 
multiple weight, otherwise, a simple one. are called the raising operators and 
the lowering ones. For an irreducible representation, there is a highest weight M, which 
is a simple weight and can be used to describe the irreducible representation. Usually, 
the irreducible representation is also called the highest weight representation and directly 
denoted by M. The Casimir C2(M) can be calculated by the formula (e.g. see (1.131) in 



N 

C2(M)=M.(M + 2p)= 5] M^rf^(A-i)^,(M, + 2), (12) 

where p is the half sum of the positive roots in the Lie algebra, is the inverse of 
the Cartan matrix, and are the half square lengths of the simple roots. 

The orbital wave function in D-dimensional space is usually expressed by the spherical 
harmonic ^^''(x) which belongs to the weight m of the highest weight representation 
(/) = (/, 0, . . . , 0). For the highest weight state, m = (/), we have 



DI 




{2i+2N-i)\ 



when D = 2N +1 



2tt^I\ 



7r^l\{l + N (13) 

1/2 

when D = 2N, 



N-l)\ 



where N^^i is the normalization factor. Its partners ^^^•'(x) can be calculated from 
Y'(f)'*(x) by lowering operators Fn{L). The Casimir for the spherical harmonic Y^\'k) is 
calculated by Eq. (12): 

L2y«(x) = C2[(/)]F«(x), C,m = l{l + D- 2). (14) 
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The spinor wave functions as well as those for the total angular momentum are different 
for D — 2N+1 and D — 2N, and will be discussed separately in the following subsections. 

A. The SO(2Ar + 1) case 

¥oY D^2N + 1 we define 

7° = (73X1, 7" = (^^2) X l<a<27V+l, (15) 

where aa is the Pauli matrix, 1 denotes the 2^-dimensional unit matrix, and (2A^+ 1) 
matrices (5a satisfy the anticommutative relations 

PaPb + (^bPa = 2(5„a, a, 6 = 1, 2, . . . , (2Ar + 1). (16) 

The dimension of matrices is 2^. Thus, the spinor operator Sab becomes a block 
matrix 

Sab^l^ Sab, Sab = -iW^- (17) 

The relation between Sab and Sab is similar to the relation between the spinor operators 
for the Dirac spinors and for the Pauli spinors. The operator k becomes 

K^a^XK, « = -i E MbLab + {D- l)/2. (18) 

a<b 

The fundamental spinor x(m) belong to the fundamental spinor representation (s) = 
(0, ... ,0, 1). From Eq. (12) the Casimir for the representation (s) is calculated to be 
C2[{s)] = {2N' + N)/4. 

The product of y'^^(x) and x(m') belong to the direct product of two representation 
(l) and (s), which is a reducible representation: 

(Ox(s)~(Z,0,...,0,l)e(Z-l,0,...,0,l). (19) 

In other words, in order to construct a wave function belonging to the representation 
(j) = (/, 0, . . . , 0, 1) there are two different ways: the combination of Yjl^ (x)x(m') and that 
of Fi^^^''(x)x(m'). They have different eigenvalues of k. Since the system is spherically 
symmetric, we only need to calculate the highest weight state for the representation (j) 
in terms of the Clebsch-Gordan coefficients 

</.|^[,0-)(x) = F(S;)(x)x[(5)] = Nn,ir-\x' + txJxM, 
\K\ = C2m-C2m-C2[{s)]+N^l + N. 

<l>-\Kmi^) = E Y^^'\±)x[{j) - m]((Z + 1), m, (.), (j) - (j)) 

m 

= ND,ir-'-\x' + ix^y {x'^'+^xiis)] + (x^^-^ + ia;2^)x[(0, . . . , 0, 1,T)] 

+ (x'N-3 + ,a;2A.-2)^[(o, . . . , 0, 1,T, 1)] + . . . (21) 

+{x^ + ^a;4)x[(l, T, 0, . . . , 0, 1)] + {x^ + zx^)x[{T, 0, . . . , 0, 1)]} , 
-\K\ = C^m - C2[(/ + 1)] - C2M +N^-l-N. 
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The wave function \l/;^_(j)(x) of the total angular momentum belonging to the irre- 
ducible representation (j) can be expressed as 

/fi(J)^K,0)(x) = /^x,0)(x), /7a.(J)^k,0)(x) = ^x,0)(x), (22) 
H^iJ)^K,u)i^) = 0, 2</i<iV-l, 
«:^K,o-)(x) = K^K,U)i^), K = ±{1 + N). 

Its partners can be calculated from it by the lowering operators F^{J). 
The radial equation will depend upon the explicit forms of (3a matrices. We express 
(3a matrices by direct products of Pauli matrices cTq 



m-1 N—m 



P2m-1 = 1 X . . . X 1 XO-i X 0-3 X . . . X O-g, 

m-1 N-m 



(23) 



/?2m = 1 X . . . X 1 XO-2 X CTg X . . . X (T3, 

P2N+1 = 0-3 X 0-3 X ... X 0-3. 
In terms of the explicit forms of (3a, we obtain 

2Af+l 

■ x) 0A',(j)(x) = r-^ 0if,O)(x) = 0-x,(i)(x), 

2N+1 (24) 
a=l 

Substituting "^fxij) (x) into the Dirac equation (5) we obtain the radial equation 

+ —G{r) = \E- V{r) - M]F{r), 
ar ^r) ^ ^p^^^ = [E - V{r) + M]G{r). 



dr r 



B. The SO{2N) case 



As is well known, the spinor representation of S0(2A^) group is reducible and can 
be reduced to two inequivalent fundamental spinor representations (+s) = (0, 0, . . . , 0, 1) 
and (— s) = (0, 0, . . . , 0, 1, 0). From Eq. (12) the Casimir for both spinor representations 
are calculated to be C2[(±s)] = (2A^^ — A^)/4. In terms of matrices given in Eq. (23), 
we define matrices for D = 2N: 

7° = /32Ar+i, r = P2N+iPa, l<a<2N. (26) 

7° is a diagonal matrix where half of the diagonal elements are equal to +1 and the 
remaining to —1. Because the spinor operator Sab and the operator k, are commutant 
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with 7°, each of them becomes a direct sum of two matrices, referring to the rows with 
the eigenvalues +1 and —1 of 7*^, respectively. The fundamental spinors x±{^) belong to 
the fundamental spinor representations (+s) and (— s), respectively, and satisfy 

7V±(m) = ±x±(m). (27) 

The product of Y^\-x.) and x±(m') belong to the direct product of two representation 
(l) and (±s), which is a reducible representation: 

(/) X i+s) ^ (/, 0, . . . , 0, 1) © (/ - 1, 0, . . . , 0, 1, 0), 
(Ox(-s)~(i,0,...,0,l,0)©(i-l,0,...,0,l). 

There are two kinds of representations for the total angular momentum: the represen- 
tation (ji) = (Z, 0, . . . , 0, 1) and the representation {j2) = 0, . . . , 0, 1, 0). Their Casimirs 
are the same: 

C2[{ji)] = C2[(j2)] = 1(1 + 27V - 1) + (27V2 _ 7V)/4. (29) 

There are two different ways to construct a wave function belonging to the represen- 
tation (ji): the combination of Yjl\-x)x+{ni.') and that of yj^+-'^^(x)x-(m'). Due to the 
spherical symmetry, we only calculate the highest weight state for the representation (j'l) 
by the Clebsch-Gordan coefficients: 

0x,(,,)(x) = F(S;)(x)x+[(+s)] = ND,ir-'ix'+tx^yx4i+s)], 

<I>-K,in)i^) = E Y^'''\^)X-[{ji) - m]((/ + 1), m, (s), (j,) - (j,)) 



m 



= ND,ir-'~\x' + ixj {(a;2^-i + ia;2^)x-[(-s)] 

+ (a;2^-3 + ix^^-^)x.[{0, . . . , 0, 1, T, 0)] (30) 

+ (a;2JV-5 + ,a;2Af-4)^_[(o, . . . , 0, 1, T, 0, 1)] + . . . 
+ {x^ + z,T'^)x-[(l J, 0, . . . , 0, 1)] + {x' + tx')x-[{^, 0, . . . , 0, 1)]} , 
K = C2[(ji)] - C2[(0] - C2[{+s)] +N-l/2^l + N- 1/2. 

For the representation (J2) = (/, 0, . . . , 0, 1, 0) we have 
<PK,U,){^) = E Y^+'\y:)x4ih) - m]((Z + 1), m, {+s), (h) - m\{j2), (^2)) 

m 

= ND,ir-^~\x^ + ixj {(x^^-i - ix^^)x+[{+s)] 
+(x2^-=^ + ix2^-2)x+[(0,...,0,l,0,T)] 

+(x2^-^ + ix'^-')x+m . . . , 0, 1, T, 1, 0)] + . . . _ (31) 
+{x^ + ix4)xH-[(l,T, 0, . . . , 0, 1, 0)] + {x^ + ix^)x+[ih 0, . . . , 0, 1, 0)]} , 
0_x,o.)(x) = y(l'/(x)x-[(-s)] = Nn,ir-\x' + ix^Yx-iis)]: 
K = C2[(j2)] - C2[(l + 1)] - C2[(+s)] + N-l/2^-l-N+ 1/2. 
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In terms of the explicit forms of Pa we obtain 

2N 



(/3 • x) 0K,o-^)(x) = r 1 ^ /3„x"0K,o.)(x) = 0_i^,(j-^)(x), 

(32) 

a=l 



a=l 

2N 



(J — 1 or 2. 



The wave function '$K,{j^){'^) of the total angular momentum belonging to the irre- 
ducible representation (j^) can be expressed as 

^lKUn)i^, t) = r-^+V2e-^^* {F(r)(/>|;,|,(,,)(x) + iG(r)0_|;,|,o,)(x)} , 
= r-^+V2e-^^t {F(r)0_|;,|,(,,)(x) + iG'(r)0|^|,(,,)(x)} , 

/ + A^-l/2, when a; = 1 



-Z-A^ + 1/2, when u^2, 



(33) 



/^i(J)*i^,0.)(x) = Z*K,o.)(x), 

i/iv-i(J)*K,Oi)(x) = 0, i/Ar(J)*if,oo(x) = *if,Oi)(x), 
i^iv-i(J)*i^,0.)(x) = *K,o,)(x), //Ar(J)*if,o,)(x) = 0, 

Its partners can be calculated from it by the lowering operators F^{J). Substituting 
\I'x(j„)(x) into the Dirac equation (5) we obtain the radial equations, which are in the 
same forms as those in D = 2N + 1 case: 

+ ^G{t) ^[E- V{t) - M]F(r), 
fF(r) (34) 
^' I + —F{r) = [E- V{r) + M]G{r). 



dr r 

IV. SOLUTIONS TO THE RADIAL EQUATION IN D+l DIMENSIONS 

Although the wave functions and the eigenvalues K are different for the D = 2N + 1 
case and the D = 2N case, the forms of the radial equations are unified 

^^^^ + -GKE{r) = [E- V{t) - M]FKE{r), 
dr r 

_dF^ + HfMr) = IE - V(r) + M]G..(r), (^5) 

dr r 
K ^±{2l + D-l)/2. 

For definiteness we discuss the attractive Coulomb potential 

y(r) = -^, C^Za>0, (36) 
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where a = 1/137 is the fine structure constant. It is easy to see that the solution for the 
repulsive potential can be obtained from that for the attractive potential by interchanging 



KE 



K-E, 



V{r) 



-V{r). 



(37) 



From the Sturm-Liouville theorem there are bound states with the energy less 
than and near M for the attractive Coulomb potential and with the energy larger than and 
near —M for the repulsive potential, if the interaction is not too strong. It is convenient 
to introduce a dimensionless variable p in Eq. (35) for bound states: 



p = 2rVM2 - E^, 0< E < M. 
Solving -F(p) from Eq. (35), 



Fke{p) 



1 M-E ^ 



2 V M + E p 
we obtain a second-order differential equation of Gke{p)- 



dGKEjp) .K 

-, 1 <^ke{P) 

dp p 



(PGke{p) 



dp" 
+ 



+ 



1 -^^ + K 



+ 



P' 



\M -E 


-1 


/ M + E 





pVM2 - ^2 



-, 1 <^ke{P) 

dp p 



Gke{p) 
0. 



From the behavior of Gke{p) at the origin and at the infinity, we define 



GKE{p)=P^e-''/'R{p), 
1 IM-E 



2i \ M + E- 



- ■ 



and obtain 



^, d^R{p) 



+ 



up^ - {2Xlj + 1) p + 2A + 1 



dR{p) 
dp 



+ [iu{X - t)p + uj{K + A) + r - A - 1/2] R{p) = 0. 



(3^ 



(39) 



(40) 



(41) 



Eq. (41) can be solved by the power series expansion method for for (3+l)-dimensions 
2T| and (2+l)-dimensions [^. The results in D + 1 dimensions are as follows 



Fke{p) 
Gke{p) 



{M±E)ET{n' + 2\ + l) , „^/2 

p e 



(M2 - ^2)1/4 

r(2A + l) ^ 2M^t{K + TM/E)n'\ ^ (42) 
X [{K + tM/E) iFi{-n',2X + l,p)Tn' iFi(l - n', 2A + 1, p)] , 



{\Fke{p)\'+\Gke{p)\') dr 



(43) 
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n 



r-A = 0,l,2,.... (44) 

where iFi(a,/3,p) is the confluent hypergeometric function. When n' = 0, i^T has to 
be positive. Introduce the principal quantum number 



n 



\K\ -{D- 3)/2 + n' = \K\ - {D - 3)/2 + r - A = 1, 2, 



(45) 



The principal quantum number n can be equal to 1 only for K = {D — l)/2 and equal 
to other positive integers for both signs of K. The energy E can be calculated from Eqs. 
(40), (44) and (45) 

1 + 



E = M 



UK' - e + n-\K\ + {D- 3)/2)2 

Expanding Eq. (46) in powers of we have 
r ^2 ^4 



(46) 



2[n+(D-3)/2]2 2[n+(D-3)/2]^ 



n+{D-3)/2 3 
\K\ 4, 



(47) 



where the first term on the right hand side is the rest energy M (c^ = 1 in our con- 
ventions), the second one coincides with the energy from the solutions to the Schrodinger 
equation, and the third one is the fine structure energy, which removes the degeneracy 
between the states of the same n. 



V. CONCLUSIONS 

In this paper we generalized the Dirac equation to (D+l)-dimensional space-time. The 
conserved angular momentum operators and their quantum numbers are discussed. The 
eigenf unctions of the total angular momentums are calculated for both odd D and even 
D cases, respectively. The unified radial equations for a spherically symmetric system are 
obtained. The radial equations with a Coulomb potential are solved by the power series 
expansion approach. The exact solutions are expressed by the confluent hypergeometric 
functions. The eigenvalues as well as their fine structure energy are also studied. Our 
solutions coincide with those in 3-1-1 dimensional |^ and 2-1-1 dimensional |2^ space-time. 
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